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^ 1 Abstract 

C - 

In this paper, we rigorously study an order 2 scheme that was previously proposed by some of the 
authors. A slight modification is proposed that enables us to prove the convergence of the scheme 
while simplifying in the same time the inner iteration. 
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1. Introduction 

In 1935 Landau and Lifschitz proposed an equation that models the magnetization in a fer- 
romagnetic material II 1311 . Supposing that the three dimensional ferromagnetic sample occupies 
^ ', some domain Q. c IR 3 and calling m the direction of the magnetization, the Landau-Lifschitz- 
^ \ Gilbert (LLG) equation reads 



O 



X 



d,m - am x d,m = -y m x H e ff in Q, 

d„m = 0onda () 



o 

C\l ■ The parameters in the equation are the damping parameter a and the gyromagnetic constant y 
The so-called effective magnetic field H eff is given by the functional derivative of the micromag 
netic (free) energy S, more precisely 



H eff (m) = -— = d 2 Am + H d (m) + H ext + Q (e • m) e (2) 
am 



where the energy £ (see [U3I.I8U1 IIP is given by 

£(m) = \ ( d 2 f |Vm| 2 dx - f H d (m) • mdx - 2 f H ext mdx-Q f (e • m) 2 dx] . (3) 

2 \ Jn Jn Jn Jn / 

The four contributions to the effective field in © and the energy in (0, respectively, correspond 
to the so-called exchange, stray-field, applied and anisotropy field or energy, respectively. The 
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material constants in © and © are the exchange constant d, the anisotropy constant Q and the 
anisotropy direction e (also called the easy axis). Furthermore, the vector field H ext models an 
applied magnetic field. We will also use the notation H an i so = Q (e • m) e. The stray field H d (m) is 
the magnetic field induced by the magnetization distribution m via the following (subset of) static 
Maxwell equations 

" curl H d (m) = in R 3 
div (Hd(m) + m) = in IR 3 . ( -* 

Below the Curie temperature, the magnetization can be described by a directional field that we 
rescale to be of unit length. It is straightforward to check that the magnitude of the magnetization 

|m(x,OI = l (5) 

is conserved by the dynamics ©. Take note that the gyromagnetic term is a conservative term 
while the damping term leads to the following energy dissipation law 

4-S(m(t)) = -— f |^m| 2 dx. (6) 



Rescaling time and redefining a allows to assume that y = 1. 

The numerical approximation of solutions to © is an important issue in applications. Nowa- 
days, numerous strategies exist in the literature - among them only few reliable ones. Classical 
schemes are based on finite differences that, as usual, are well adapted to Cartesian grids. On the 
other hand, finite elements approximations are well suited in case of complex geometries and weak 
solutions, though bearing the drawback that they are in practice difficult to analyze. In particular, 
proving the convergence of a finite element solution towards a solution of © as the space and time 
steps tend to zero turns out to be quite difficult and has probably been first established in [4]. This 
result was further improved in [0] and fly, for the case where only the exchange term is present. 
We hereafter study a further generalization of the scheme proposed in yj]: An order 2 (in time) 
variant. Numerical tests support the performance of the method. 

Let us start with brief outline of our paper. In Section |2] we first recall the notion of weak 
solutions. Section [3] introduces the finite elements spaces. Section @]restates the order one scheme 
as proposed in [Q]]. The nonlinearity of the LLG equation calls for recurrent renormalization of the 
time-discrete approximation. This issue is also discussed in Section [51 Section |6] finally provides 
a derivation of our new scheme, the main result about its convergence and its proof. 

2. Notion of weak solutions to LLG 



Let us recall the notion of a weak solution to © from [5] and IU6I1 



Definition 1. Consider an initial magnetization, i.e., a vector field m e H l (Q) 3 that is a.e. of unit 
length. A vector field m is called a weak solution to © with initial data m if for all times T > 
there holds 

1. me H l (Q. T ) 3 with Q. T = Q. x (0, T), and |m| = 1 a.e. 

2 



2. for all test functions *P 6 H l (Q. T ) 3 



= d 2 Yj \ ( m x ^, m ) • dxF dx dt - I m x (H rf (m) + H CT( + H anho (m)) -W dxdt, 



^•Tdxdr-tt I (m x d,m) ■ T dx dr 



(7) 



3. ?/ie magnetization initially satisfies m(x, 0) = mo(x) m the trace sense, and 

4. the energy decreases according to 



a f \d, 

JQ.T 



S(m(T)) + a |d r m| 2 dx df < £(m(0)). (8) 



3. The finite element scheme 



As in [4], our discretization relies on piecewise linear finite elements in space combined with a 
linear interpolation in time. The domain Q. is discretized by a conformal triangulation T h of mesh 
size h with vertices (x' ! )i<,<w,,. Let us denote by (0- 1 )i<,<^ the set of associated piecewise linear 
basis functions that satisfy 0- ! (x'O = 6(j at the vertices x h . for 1 < i,j < N/,, where Sy denotes the 
Kronecker symbol. This amounts to a standard i ,1 (7 _ /,)-discretization. Based on the scalar basis 
{<p h i)\<i<N h we construct the vector-valued finite element space in the form of 



Due to the constraint ©, the solution to © is sought for in the subset 

M h = [u h G V h , s.t. Sfi, u, e § 2 } c V h . 



Let us also introduce the tangent space in m h = m;0; £ is denoted by 



K m = | v/! = Z S - t W ' V >' ' m < = °J 



Furthermore, the classical nodal interpolation operator is given by 

J A :C°(ara 3 )-> V h 

u^^u(xf)^. (9) 

To simplify notations, the index h of the ansatz functions will be neglected from now on most of 
the times, i.e., we write u, v, etc. instead of u\ v\ respectively, in case this does not lead to any 
ambiguities. 

3 



4. Revisiting the 0-scheme 



The finite element scheme proposed in |4J] relies on the observation that the LLG equation (QQ) 
- with the notation v = d t m - can be rewritten in the following weak form 



a 



J v • ¥ dx + a J 

Jn Jn 



m x v *P dx 



= -d l 



Jn 



Vm VFdx + 



f (H 

Jn 



d (m) + H ext + H amso (m)) T dx. (10) 



Equation (PTOl) holds for every test function *P e H l (Q,, R 3 ) that satisfies *F(x) • m(x) = for a.e. x 
in Q.. The reformulation of (OQ) in the form of (TT0T > motivated the following 6- scheme introduced 
in [Si: 



Algorithm 1. Given an initial m e M/, choose 6 e [0, 1] an J a rime ste/? size r = jj wjY/z iVeN. 
Forn = 0,1,..., iV 

a J ^znJ v" 6 K m n such that for all test functions *F e 

a ( v"-Tdx+ ( m"xv"-Tdx 
Jn Jn 

= -d 2 \ V(m" + 0rv") • VP dx + | (H d (m") + H eAt + H aniw (m")) T dx < 1 1 > 

Jn Jn 



m" + rv" I 



It is noteworthy that this procedure requires the solution of a linear equation in each time step only. 
Moreover, due to the fact that the symmetric part of the underlying matrix is positive definite, 
existence and uniqueness of a solution to (TTTT) is guaranteed. 

rr 

The time discrete solution constructed via algorithm (TTTb at time-steps N - 
as follows: 



is interpolated 



Definition 2. In each time interval t e [nr, (n + 1 )r) with n e {0, • 

t - nr 



, Af} we set 



ni/, r = 



n+1 , (/z + 1)t - f „ 
m h m 



T 

m", v /liT 



Our notational convention is thus that ni/ ) T , m^ r and V/, jT refer to suitable time interpolants of the 
time discrete approximation m" and v". Notice that m^ jT is piecewise linear in time whereas mj ir 
and V/j >T are piecewise constant. (The introduction of the piecewise constant magnetization will 
be useful in the convergence proof.) Based on this discretization, weak convergence of the con- 
structed approximation was established in QJLD - Both the proof of this result and the proof in case 
of our new scheme consist of the following two main "classical" steps: As a first step establishing 
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an energy estimate which guarantees the convergence (sufficiently strong) of the sequence con- 
structed and then in a second step verifying that the limit indeed satisfies the equation. As far as 
the first step is concerned, the following section addresses the fact that the energy behaves well 
under renormalization - in principle a strongly nonlinear modification of the flow. 



5. Renormalization decreases the energy 



The influence of the renormalization on the exchange energy was for instance investigated 
in [0] in the continuous case. More precisely, it was shown that for maps w G H l (Q, M 3 ) with 
|w(x)| > 1 a.e. x e Q. one has 



f V^ 2 dx< f 
Jn Iw| Ja 



dx < |Vw| 2 dx. 



(12) 



Hence, the renormalization step is expected to be energy decreasing - a least as far as the Dirichlet 
energy is concerned. Applications more related to finite element approximation of micromagnetic 
configurations can be found in yfl. The discrete version of (fT2l was proved by Bartels in D6fl: 



Theorem 1. [6] If the basis functions of the P l -approximation satisfy 

f V0f-V05dx<O, 
Jn 

then for all v = £i v,-^ 1 G Vh such that Vz G { 1 , • • • , Nh), |v,-| > 1 it holds that 

J 

Jn 



(13) 









dx< r 




Jn 



(14) 



In 3d , the condition (PT31) - and hence (PT41) - is for instance satisfied provided all dihedral 
angles of the tetrahedra of the mesh are smaller than n/2, see iflill . 



6. The new (almost) order 2-scheme 

Let us embark on the motivation and description of our new scheme. As remarked in lU], it 
is not sufficient to choose 6 = \ in (fTTT) to achieve quadratic order due to the renormalization 
which inherently introduces an error of order 2. Hence, it is necessary to modify the time-discrete 
approximation of the magnetization m. 

Consider an iterate m(nr) at time nr. It is well known that the mid-point rule is exact up to 
cubic error, i.e., 

m((n + 1) t) = m(nr) + nn,((n + \) r) + 0(r 3 ). 
Now, given a current iterate m(nr) at time nr, a Taylor expansion up to cubic order, i.e., 

T 2 

m((n + 1) t) = m(nr) + Tm ? (nr) + — m tt (nT) + 0(t ) 

5 



reveals that the parallel component of the subsequent iterate (along m(nr)) is due to the unit length 
constraint given by 

m((n + 1) t) • m(nr) = 1 - T 2 \m t (nT)\ 2 + 0(r 3 ). 
This can easily be inferred from the unit length constraint by differentiation, i.e., using the relations 

m m, = 0, 
m, m + m m = 0. 

We therefore propose to modify the original first order scheme by replacing the tangential update 
with the following higher order approximation 

v = P m ± m,((w + i) t) 
= P m ± (m,(nr) + |m tt (/rr)) + 0(t 2 ) 

= m ; (nT) + \ P m xm„(nr) + 0(r 2 ), (15) 

where i> m ± denotes the projection onto the orthogonal component of m(nr). 

We will use the short hand notation m = m(nr) and m, = m t (nr) - provided that what is stated 
remains clear without ambiguity. Let us proceed with the derivation of the equation that is satisfied 
by v = m f (nr) + 5 P m ±m n (nT), i.e. the counterpart to <TT0T > . The equation will be inferred from the 
differentiated LLG equation which we restate as 

am, + m x m = H eff (m) - (H eff (m) -111)111 (16) 

by multiplying (Q~|) with mx. 

Remark 1. Although the mid-point rule is of order 2, our scheme will be only almost of order 2 - 
as the section 's title suggests and as we will see in the sequel. We have to introduce a regularizing 
term in order to obtain the necessary estimates in the convergence proof. This term prevents the 
scheme from being of order 2, in the sense that the consistency error is not of order 0(t 3 ) but only 
0(r 3_e )/or any e > 0. On the other hand, this regularization approach allows for unconditional 
convergence of the scheme. If we do not insist on unconditional convergence, then under the 
condition r «.h, consistency up to order 0(r 3 ) is attainable. 



To begin with, the differentiation of (1161) w.r.t. time yields 
am„ + m t x m„ (17) 
( IIeft (m,) - |^^(m t ) m) m - (H eff (m) m ) m (H eff (m) • m) nr . (18) 



dm \ dm 

where 



= d 2 Am, + H d (m,) + Q (e • m f ) e 

dm 
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and where we once again used the unit length constraint ©. The application of the projection to 
(PT8l) in combination with (fT6l) yields 

ay ■ *P + m x v *P dx 

= £ H eff (m) • T dx + \ £ ^(m r ) ¥ dx - \ J"(H eff (m) • m) m, ■ *F dx 

for any test function *P with ¥ ■ m = 0. Observe that m,(nr) = v + 0(t), cf. (fT5l) . Therefore up to 
higher order terms 



Jn 



f («+|(H, 
Jn 



eff(m) • m)) v 'P + mx 



v Ydx-f f 

Jn 



eft- 



^m 



(v) -^dx 



Jn 



H eff (m)-^dx + 0(T 2 ), (19) 



where we remind that m = m(nr) and m t = m ( (nr). Observe that the latter equation is (at first 
sight surprisingly) linear in v. However, nothing can be stated about its well-posedness since 
both the first and the last contribution on the l.h.s. of (fT9l potentially affect the definiteness of the 
symmetric part of the operator. In order to guarantee solvability and uniqueness we proceed with 
higher order modifications that will finally lead to a well posed formulation. We address the first 
contribution and define 



J a + \ min(x, M) for x > 0, 
<Pm(x) = i a 



1 + f min(-x, M) 



for x < 0. 



Notice that <Pm(x) = ot + \ min(jc, M) + 0(t M ). By abuse of notation we define 

(p M (m) = <£ M (H eff (m) • m). 



(20) 



(21) 




Figure 1: The regularizing cut-off function ip(x). 

As long as H e ff(m) • m is uniformly bounded, we derive from (TT9T > by plugging in (|2TT> that 

f cp M (m) v • Y + m x v Ydx - f f ^P^(v) Ydx = f H eff (m) • Tdx + <9(r 2 ). (22) 
Jn Jn dm J n 



Replacing H e ff and by their very definition, we obtain the counterpart to (fT0T > for our new 



second order scheme: 
^m(hi) v • *P + m x 



Jn 



Jn 

Jn 



d 2 VvVP- H d (v) • *P - <2(e • v)(e • *P) dx 



- d 2 Vm • V*P + H d (m) • *P + Q(e ■ m)(e • *P) + H ext • ¥ dx. (23) 



We introduce only one further, final modification which implements the strategy delineated in 
Remark \T\ In order to maintain unconditional convergence we additionally modify the second 
highest order term on the r.h.s. in the following way 



if 

Jn 



d VvVPdi 



if 

Jn 



(1 +p(r))J 2 Vv VTdx, 



where p(r) — > as r — > 0. Take note that for p decreasing at least linearly, quadratic order is 
conserved. However, only in case that p is slightly sublinear, for example p(r) = r| ln(r)|, do we 
in fact achieve unconditional convergence. 

Adopting AlgorithmCD we arrive at the following scheme: 

Algorithm 2. Given an initial m° 6 M h choose a time step size r = ^ with N e N and appropriate 
p(r) and M, cf. Theorem\2\ For n = 0,1, . . . ,N 

a) find v" 6 such that for all test functions T e K m « 



L 

Jn 



<p M (m n ) v" • T + m" x v" • *F dx 



+\ \ (1 + p(r)) d 2 Vv" • VP - H d (x n ) • T - <2(e • v")(e • dx 

_ Jn 

- d 2 Vm" • + (H,(m") + H ext + H ami0 (m")) • Tdx. 

/n 



X 



(24) 



b)setm" +l = ^]m;. !+ Vf, where \/i, m^ +1 



m" + rv" 
|m" + rv" I 



The appropriate choice of p and M can be inferred from our convergence result, see Theorem 

Let us sum up: The new scheme replaces the search of v as solution to (fTTT) by the search of 
v as a solution to (1241) . Besides this substitution, the algorithm outlined in Section |4] remains as 
before in the sense that the renormalization and the interpolation w.r.t. time are left unchanged. 
Since equation (1241) is linear in v, our algorithm is very favorable in practice. 

Before we state our theorem about the convergence let us explicitly make a statement about its 
order. 



Proposition 1. Consider a smooth (in space and time) solution m to (l2"4l at time t + r and a semi- 
discrete time-approximation to m at time t + t on the basis of (1241) . More precisely, given m at 

8 



time t = nr determine v = m,(nr) + | P m im tt (nr) as a solution to the variational formulation (1241) 
wzY/j p(r) = and M(t) sufficiently large and set 

m(x, t) + tv(x, 

m(x, ? + r) = for all x e Ll. 

|m(x,0 + rv(x,0r 

TTzen m(? + r) approximates m(t + r) ap to cwozc error in r. 

Argument for Proposition^ The proof is a direct consequence of the Taylor expansion performed 
in© 



Remark 2. 77ze smoothness of solutions to © /las oeen widely studied during the course of the 
past years. In general, the formation of singularities cannot be ruled out and we can usually not 
assume that a solution to the initial value problem will be regular. Our statement about the order 
of the approximation is thus only a first little step on the way to a proof of the order of convergence, 
which is way beyond the scope of this paper. 

Let us now turn to the convergence result. 

Theorem 2. Let m e H l (Q,,S 2 ). Suppose m° —* mo in H l (£l) as h — > 0. If the regular sequence 
of conformal triangulations (7~h)h>o satisfies condition (fT3l) . then the approximation (m^ jT ) of the 
sequence constructed via Algorithm\2\and interpolated according to Definition\l\converges (up to 
the extraction of a subsequence) weakly in H l (£l T ) to a weak solution m of (OQ) as h and r tend to 
provided p(r) — » T _»o and one of the two following conditions hold: 

• t"V(t) -»(a, t )->o 00 and tM -> (ft , T) ^o or 



p = and t <sc h as (h, r) — > 0. 



Proof of Theorem\2\ As stated before, the proof consists of two main steps: Establishing estimates 
which guarantee the existence of a sufficiently strong converging subsequence, and finally proving 
that the latter converges indeed to a solution (which satisfies the energy estimate). We will need 
the following classical estimate from elliptic regularity theory, namely 

l|H d (m)|| iP(n) < Cllmll^n), (25) 

for all p 6 (1, +oo) and for positive constant C which depend only on p. 

Bounds on the sequence. As we have already observed, the variational formulation in the iteration 
of (1231) possesses a unique solution v". We test the equation with *P = v" itself to find that 

f (f M (m n ) |v'f dx + \ f(l + p(r)) d 2 |Vv'f - H d (v") • v" - Q(\ n ■ e) 2 dx 

- d 2 Vm" • Vv" + H d (m") • v" + Q(e • m")(e • v") + H ext • v" dx. (26) 



I 

Jn 
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Since we assume that the triangulation TJ, satisfies the angle condition (1141) we have that 

f |Vm" +1 | 2 dx < J |V(m" + tv")| 2 dx 
Jn Jn 

< f |Vm M | 2 dx + 2t f Vm" • Vv" dx + t 2 f |Vv"| 2 dx. 
Jn Jn Jn 

Using (1261) we obtain that 

J 2 f |Vm" +1 | 2 dx < d 2 f |Vm"| 2 dx-2r f ^(m") |v'f dx + r 2 f H d (v") • v" + Q(e • v„) 2 dx 
Jn Jn Jn Jn 

+ 2r f H d (m")-v' l + H aniso (v' 1 )-v' 1 +H ext -v"dx-T 2 p(T)J 2 f |Vvf dx. (27) 
Jn Jn 

Before we move on, let us just rewrite the latter estimate as 

|Vm M+1 | 2 dx< d 2 f |Vm n | 2 dx-2r f w M (m n ) |v'f dx + r 2 f — ^(v") • v" dx 

Jn Jn dm 

(28) 



n ^n 



+ 2r f H eff (m n ) • v" dx - r 2 p(r) J 2 f |Vv"| 2 dx 
Jn Jn 



We partially neglect the negative contributions on the r.h.s. of (12~ST) - those which are quadratic in 
v" - and use (1251) to obtain 



(I 2 f |Vm" +1 | 2 dx 
Jn 



< d 2 f |Vm'f dx - 2r f ^ M (m) |v"| 2 dx + 2r||H eff (m' ! )|| L 2 (f2) ||v"|b (n) - r 2 p(r) J 2 f |Vv"| 2 dx 

Jn Jn Jn 

< d 2 f |Vm'fdx-2r f ^(m) |v"| 2 dx + Cr|| Vn || L2(n) - r 2 p(r) J 2 f |Vv"| 2 dx, (29) 

Jn Jn Jn 

where the generic constant C depends on Q and |Q|. Due to Young's inequality, we have that 



Cr||v„|| L 2 (n) < r^||v n ||2 2(n) + for/? > 0. Using the uniform bound 

we find by rewriting (|29l that 

d 2 J |Vm M+1 | 2 dx +^r||v"|| 2 2(n) + r 2 p(r)d 2 J |Vv"| 2 dx 

<^ 2 fW dx + TC2(g ' |Q|) . (30) 
Jn 4/3 



10 



Summing up in (|30l) over the time steps we find that 



d 1 | \Vm N \ 2 dx+firV f |v"| 2 dx + T 2 p(r)d 2 f |Vv'f dx 



.ji|2 



?rV f IVf 



I 



< C T, J 2 |Vm | 2 dx,j8,e,H a 



(3D 



From now on, most of the arguments follow the same line as in QlJ]. It holds that 

< |v"|, for all n < N, audi e {!,■■■ ,N h }. 



m" +1 - m 



Moreover, there exists c > such that for all 1 < p < +oo and all (p h e Vj, there holds 

c 



<^^|0,(4)| p <c|| 



LP (SI)' 



which implies 



m" +l - m 



< c 2 ||v"|| i2 



L 2 



Hence we obtain from the energy estimate (13TT ) using d33T > the following bounds 

m ft T is uniformly bounded in H Y (Q, T ), 
\h,r is uniformly bounded in L 2 (Q. T ). 



(32) 



(33) 



(34) 
(35) 



Due to (1341) and (1351) . there exist m e H l (Q. T ) and v e L 2 (Q. T ) such that up to the extraction of 
subsequences 



mh.T -^(fc,T)->o m weakly in H\Q T ), 
m/,, r -»(ft, T )-»o m strongly in L 2 (£2 r ), 
v/,,r -^(A,r)^o v weakly in L 2 (Q T ). 

In addition, we have from (13TT) that 



(36) 
(37) 
(38) 



V r 2 p(r) f IVv'f dx = rp(r) f f |Vv,, T | 2 dx < C 
^ Jn Jo Jq 



< +oo 



If p decreases only sublinearly, i.e. r ! p( T ) ~~ *t->o we deduce that 

T||Vv|| L 2 (nr) ->(/,, T )-»0 0. 



(39) 



If p decreases linearly or faster we have to resort to the inverse estimate ||Vv|| L 2 (nr) < ^||v]] L 2 (f2r) 
in order that estimate (1391 holds true. In fact, is easily seen that (1391 is follows from the inverse 
estimate in case of t <sc /z. 
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Preliminary estimates. We want to prove that m satisfies © and follow the strategy of [U]. To 
begin with, we restate some further estimates from [Q]] and derive some necessary statements 
about convergence. Observe that for all n = 0, • • • ,J and all t e [tit, (n + l)r) 



\m hiT (x,t)-m h (x,t)\ = 



(t - m) 



m" +1 (x) -m"(x) 



<T|d f m A>T (x,?)| . 



Therefore 
which entails that 



m h T - m h T \\ L 2 iiiT) < T\\d t m htT \\ L2(nT) -»0, T )_»o 0, 



m h,r -*(h,r)^o m strongly in L 
Moreover, on any tetrahedron K of 7"/,, and for any u e M/, one has, x 1 . being any vertex of K, 

||u(x)| - |u(x?)|| 2 < C/i 2 |Vu| 2 , 
(recall that Vu is constant on K), from which one deduces (since |m^ r (x' 7 )| = 1) 

f |l-K T l| 2 dx<0 2 ||Vm,- r || 2 2(nr) . 

Jci T 

This shows that |m(x, 01 = 1 a - e - ( x > t) e Q. T . 

Eventually, from the fact that at each vertex Vi e {1, • • • , Nf,} 



Imf 1 - m" - rv" 



K + r\'J\ - 1 < {T 2 \y1\ 2 , 



(40) 



we derive 



m M+1 - m" 



Appealing to (1321) the latter entails that 



— 2 1 ( 1 



9 t m KT - V A>T || il(0r) < C T"||V/ !iT || i2(ar) -^(h,r)^0 0. 

This is sufficient to conclude that v = d t m in (1381) . 

General properties of interpolation operator. Before we start with the penultimate step of proving 
convergence, let us state some general properties of the nodal interpolation operator which we 
repeatedly use in the sequel. Up to dimension three, there holds for any function tp e H 2 (£l) c 
C°(Q) 



\\(p - IhiMnKm < Cfc||VVlbn- 
Since the basis functions are linear on each triangle one can deduce form (I4TT) that 

seeflp.7]. 
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(41) 



(42) 



Convergence to a solution of the LLG equation. Having established the preliminary results above, 
we are now ready to proceed with the proof of convergence: Test (1231) with T = Ja(mj~ x *F) 
where *F e C^(Q. T ) 3 - We recall that I h is the nodal interpolation, cf. ©. After suitable integration 
in time we hence obtain from (1241 with the choice of Y = I/ 1 (m hT x that 

I (p M (mJ 1T )vi hT - I h (ml T x } ¥)dxdt+ I m£ T x v ft , T • I h (m h T x *F)dxd? 



+ 



(1 + p(T)) d 2 Vv h>T ■ V J ft (m- T xT)- H d ( V/! , T ) • J k (m£ T x *) 
- Q(e • v k>T )(e • I h (vn KT x *)) dxd? 
- d 2 Vm A T • VJ A (m£ r x *) + H d (m, ;T ) • I h (m h T x T) 

+ Q(e • m^Xe • X h (mZ T x T)) + tf ext • J k (m^ T x *F) dx dr. (43) 



Our goal is to pass to the limit (t, h) — » in the latter equation (|43l to recover the LLG equation 
(PTOl) . As we shall see, the first and the third term on the l.h.s. and the first term on the r.h.s. are a 
little bit subtle and have to be treated with caution. The remaining contributions behave well under 
the established convergence; this is particularly due to the fact that H d is L 2 -continuous. For the 
second contribution on the l.h.s. one further uses that the L°° bound on m improves (1371) to strong 
convergence in any LP with 1 < p < +oo. 

Let's start with the first contribution on the l.h.s. Observe that \(p M \ is uniformly bounded. 
Moreover it holds that \tp M - a\ < As long as tM — > for (h, r) — > the strong convergence 
of mr is sufficient to conclude that 



h,r 



I (f M (m h T ) \ Kt ■ I h (m h T x T) dx dt -»(*, T )_»o a \ 

JQ. T j£i T 

In fact, using the triangle inequality we find that 

(f M (ml T ) \ n hT ■ J A (mj~ T xV)dxdt-a I v • (m x T) dx dr 



v • (m x T) dx dt. 



(44) 



< 



<p M (m h T ) v k>T • (m h T x T) dx dt - a 



+ 



J 



(mxY) dxdt 



<p M (m h T ) y Kr ■ ((m h T x T) - I h (m h T x ¥)) dx df 



(45) 



The first term tends to zero since 



in L°°(Q), 
L 2 (^ r ), and 



m A,r -*(h,r)-»o m in 

va.t ->(*,t)-»o y = ~^f m L (Qt-) 
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as h, t — > 0. 

Since (p M (m hT ) is uniformly bounded, we can evoke (1421 to obtain that the second contribution 
tends to zero. This establishes (1441) . 

Let's turn to the next term in (I43T ). Convergence in this case essentially relies upon the estimate 
(1391) . In fact, appealing once again to (1421) we see that instead of establishing 



is f 

Jn, 



(1 + p(T))Vv ;j , r • VJ A (m^ T x dx -> ( a, T )_* 0, 



if suffices to establish 



\d 2 f (l+p(T))Vv*, T -V(m^ T x*)dx-» ( A, T) _oO, 
Jn r 



(46) 



(47) 



which follows obviously from (1391 using Young's inequality. 

Finally, the convergence of the last term in (|45T ) follows from the orthogonality property of the 
cross product and (136T >. (137T > by once again appealing to ([421 since 



Vm A >T • V J/,(m/T, T x *F) dx dt - Vm mxVTdx d? 

Jn r ' 

f V m/ ; r • V (l h (m- hT x *) - (m^ T x *)) dx d* 

Jn T 

Vm^ iT • V(mj~ T xf) - Vm • m x VP dx d? 



^3 



(48) 



f Vm,, ■ V (l h (m- hT x *) - (m,- T x *)) dx < 

Jn r 

I Vm ft r • (m£ T x VT) - Vm (m x V*) dx df 

Jn r 



(49) 



(50) 



Energy estimate. We finally establish the energy estimate. From (127T > we deduce that Vn e 
{<>,•••,#*} 



L 

Jn 



fi(m' ,+1 )-£(m") < -2ar I ^m") |v"| 2 dx + 2r I H eff (m") • v" dx 



n 



+ r 



r ^£ (v «). v « dx _ T 2 P(T)J 2 

Jn 



/n u m 



<9„ 



r ivv"r 

Jn 



dx 



f H eff (m" +1 + m") • (m' ,+1 - m") dx, 



(51) 



cf. ©. Let us introduce another short-hand notation for the remaining effective field, namely 
H" ff = H e ff(m"). We consider the contributions in f|5TT > separately and start with the observation 
that 



2r 



f H" ff ■ v" dx - f (H^ 1 + H' e ! ff ) • (m" +1 - m") dx 
Jn Jn 



= 2 f YT eS ■ (m n+1 - m n - tv") dx + f (H^ 1 - H" 
Jn Jn 
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eff )-(m' l+1 -m")dx. 



Hence due to 03]) and (gO]) combined with ((321) 



2r 



f H' e ! ff ■ v" dx - f (H^ 1 + H^ ff ) • (m" +1 - m") dx 
Jn Jn 



dx < Cr 2 (||v"|| L 2||v"|| L 4 + ||v"|| 2 2 ) (52) 



In order to bound the stray-field contribution we have employed (1231) with p = 4. The contributions 
in the second line of the r.h.s. of (15TI) are of higher order in r. The first term can be easily bounded 
using Young's inequality: 



Jn 



eft 



(v") • v" dx 



< Q\y 



n\\2 



in dm 

Plugging in ([52]) and ([53]) into (ED yields that 

fi(m" +1 ) - £(m") + 2r \ <p{m n ) |v'f dx 

Jn 

< Cr 2 (||v"|| L2 + ||v"|| 2 2 + ||v"|| L2 ||v"|| L 4 +p(r)d 2 ||Vv%) 

< C'r 2 (||v"|| L2 + Hv"!! 2 , + ||v"|| L2 ||Vv' ! |b +p(r)||Vv"||2 2 ), 

where C denotes a generic constant. Here we made use of the classical Sobolev embedding 

||v' ! || L 4 < C||Vv' ! || L2 . 

Summing from n = Oto N - \ leads to 



(53) 



I 



8(m(NT)) - S(m(0)) + | ^(m^ T )|v„, r | z dxdt 

< Ct(||v a>t || L 2 + ||v ;j , T || 2 2 + ||v ftjT |b||Vvft ;T |b +p(r)||Vv/ 1 , r || 2 2 ). 

We are now ready to pass to the limit. Noticing once again that T||Vv' ! || L 2 (f2r) ) is uniformly bounded 
from (1391 ) we derive that 



£(m(iVT)) - £(m(0)) + a 



Jo Jn 



dx dt < 0. 



(54) 
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